Revealing the nature of high-mass X-ray binaries 
through multi-wavelength and statistical analyses 
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Abstract 

We summarize the results of our long-running campaign to help understand the nature of high- 
mass X-ray binaries (HMXBs), emphasizing recent Suzaku observations of IGR J16207— 5129 and 
IGR J17391— 3021. Thanks to the expanding ranks of HMXBs in our Galaxy, we are able to perform 
more reliable statistical analyses on the three currently- known sub-classes of HMXB: those with super- 
giant companions (SGXBs); those with Be companions (BEXBs); and the enigmatic Supergiant Fast X-ray 
Transients (SFXTs). We discuss new diagnostic tools, akin to the "Corbet diagram," in which HMXBs 
tend to segregate based on their dominant accretion mechanism. We show how SFXTs span across the 
divided populations of BEXBs and SGXBs, bolstering the intriguing possibility that some SFXTs repre- 
sent an evolutionary link. The use of HMXBs as tracers of recent massive star formation is revisited as we 
present the first ever spatial correlation function for HMXBs and OB star-forming complexes. Our results 
indicate that at distances less than a few kpc from a given HMXB, it is more likely to have neighbors 
that are known massive-star forming regions as opposed to objects drawn from random distributions. The 
characteristic scale of the correlation function holds valuable clues to HMXB evolutionary timescales. 
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1. Multi-Wavelength Analyses 

High-mass X-ray binaries (HMXBs) feature a compact 
object such as a neutron star (NS) or a black hole candi- 
date (BHC) accreting material from a young and massive 
donor star (M > 10 Mq). A majority of HMXBs host 
main-sequence Be stars that have not filled their Roche 
lobe. These systems (BEXBs) are usually transient with 
flares produced whenever the sometimes wide and eccen- 
tric orbit brings the compact object close to its compan- 
ion. Supergiant HMXBs (SGXBs) are accompanied by 
an evolved supergiant O or B star whose wind steadily 
feeds the compact object. Their emission tends to be per- 
sistent with low-level variability stemming from inhomo- 
geneities in the wind. Supergiant Fast X-ray transients 
(SFXTs) are an emerging class of HMXB. They share 
properties from both established groups: like BEXBs, 
they are characterized by outbursts of short (hours) to 
long (days) duration whereby the peak intensity is sev- 
eral orders of magnitude greater than during quiescence; 
and like SGXBs, they are systems in which a compact 
object (usually a NS) is paired with a supergiant OB 
star. 

While INTEGRAL has proven adept at finding new 



HMXBs^ , in particular the heavily-absorbed or fast tran- 
sient sources, its spatial resolution is of the order of an 
arcminute, i.e., too wide to permit a search for opti- 
cal/IR counterparts in the busy fields in which HMXBs 
are located. It is only through dedicated follow-up obser- 
vations with focusing X-ray telescopes such as Chandra, 
Swift-XKT, and XMM-Newton, that an arc-second posi- 
tion can be obtained that is of use for optical/IR studies. 
Following up an IGR discovery with RXTE and Suzaku 
can provide more precise timing (RXTE), broader spec- 
tral coverage (both), and a lower background (Suzaku). 

A number of research groups have put in place cam- 
paigns to systematically follow up new or previously- 
known HMXBs (and other source types): e.g., Chaty 
et al. (2008); Landi et al. (2010); Masetti et al. (2008- 
2010). Our group has been active in the analysis of aU 
of the aforementioned telescopes: e.g., Bodaghee et al. 
(2010, 2011a); Rodriguez et al. (2009, 2010); Tomsick 
et al. (2008, 2009, 2011). Recently, we analyzed Suzaku 
observations of the SGXB IGR J16207-5129 and the 
SFXT IGR J17391-3021 (=XTE J1739-302). The ob- 
servation of IGR J16207— 5129 concludes with 30 ks of 

*^ An up-to-date list of INTEGRAL sources and their properties 
can be found at: http://irfu.cea.fr/Sap/IGR-Sources 
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Fig. 1. Reported X-ray column densities for HIVlXBs as a function of 
orbital period are presented for SGXBs (blue disks, of which the 
SFXTs among them are boxed), BEXBs (green filled squares), 
and HMXBs with unclassified companions (red filled triangles). 



emission consistent with no flux. There are no signifi- 
cant spectral changes, and the time scale is an order of 
magnitude too long for an off-state. We tentatively at- 
tribute this event to an eclipse of the X-ray emitter by 
its supergiant companion, but note that the demarcation 
between SFXTs and SGXBs is no longer as clear as was 
once believed, and so whatever mechanism is responsible 
for the quiescent state of SFXTs could be at play here. 
For IGR J17391— 3021, we observed quiescent emission 
interrupted by 3 weak flares in which the peak luminos- 
ity is only a factor of 5 times that of the quiescent phase. 
In this case, the change in luminosity was accompanied 
by a significant increase in the absorbing column indicat- 
ing that the accretion of obscuring clumps is responsible 
for the weak flares. When placed in context of the long- 
term monitoring by Swift (Romano et al. 2009), we find 
that these low-activity events (including those seen in 
this source by Bozzo et al. (2010) using XMM-Newton) 
represent the most common state for this source (60 ±5% 
of all observations in the 0.5-10 keV band). 
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Fig. 2. Reported X-ray measured column densities for HMXBs as a 
function of pulsation period. Blue disks represent SGXBs (SFXTs 
are boxed), while green squares represent BEXBs. Unclassified 
HMXBs are denoted by red triangles. In the bottom panel, we 
focus on SGXBs and SFXTs and show equilibrium spin periods for 
a wind-accreting neutron star model from Waters & van Kerkwijk 
(1989) assuming a variety of magnetic field strengths and binary 
separations. 



2. Absorption Modulation 

The X-ray photoelectric absorption (iVn) has been mea- 
sured for a number of HMXBs. The intrinsic absorbing 
column is a direct measure of the amount of material in 
the vicinity of the compact object. As one would ex- 
pect, the longer the binary period (i.e., the greater the 
average separation), the lower the iVn, leading to a weak 
anti-correlation (Spearman rank correlation i?s = —0.3, 



or a 2% probability of being due to chance). Overall, 
SGXBs and BEXBs are segregated into distinct regions 
of the plot while the SFXTs bridge the divide. The ver- 
tical scatter in the plot is noticeably reduced when the 
column densities are normalized by the interstellar value 
expected along the line of sight (e.g.. Dickey et al. 1990, 
Kalberla et al. 2005). 
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Fig. 3. The spatial distribution of HIVIXBs in Galactic coordinates with symbol size proportional to the X-ray measured column density as 
reported in the literature. The contours denote cumulative line-of-sight absorption levels of 10^^, 5 X 10^^, and 10^^ cm~^ (Dickey et al. 
1990). A cluster of heavily-obscured sources can be seen in the direction of the Norma Arm with no equivalent overdensity towards the 
other arms. 



Whether absorbing material is accreted or spun away, 
the pulsar spin frequency will be affected. When spin 
periods of the neutron stars in HMXBs are plotted as 
a function of the iVu, we observe a segregation between 
BEXBs and SGXBs. The exceptions are SGXBs such as 
Gen X-3 that have atypical accretion regimes for their 
class. The SFXTs straddle the disparate groups. A 
weak correlation is detected with i?s = +0.3 which cor- 
responds to a 1% probability of being due to chance. 
One explanation for the correlation is that the accre- 
tion of high-density wind material around the neutron 
star (with no preferred direction of angular momentum) 
will tend to slow the pulsar spin. Another explanation 
is that the perceived correlation is simply the result of 
two populations occupying different regions of the dia- 
gram. Once again, the vertical scatter can be diminished 
by normalizing the column density of each source with 
respect to the absorption expected in that direction. 

These figures lend support to the emerging model 
in which some SFXTs are descended from BEXBs and 
could therefore represent an evolutionary link between 
BEXBs and SGXBs (e.g., Negueruela et al. 2008, Ghaty 
et al. 2011). First appearing in Bodaghee et al. (2007), 
these diagrams have been updated to distinguish the 
SFXT class, and to include the latest results such as 
MAX! J1409-619, an unclassified HMXB with a spin 
period of '-~^500s and A^h ^ 3 x 10^^ cm^^ (Kennea et al. 
2010ab) 



Figure 2 (right panel) shows the spin period vs. A'h 
diagram with lines of equilibrium spin period as derived 
from the wind-accreting pulsar model of Waters et al. 
(1989). Ganonical values are assumed for the mass of 
the neutron star and its supergiant companion, as well 
as for the supergiant wind velocity and the neutron star 
accretion rate. The column density for the model is ob- 
tained by integrating a range of wind densities along the 
line of sight. Several binary separation values and neu- 
tron star magnetic field strengths are considered. Pulsars 
in these HMXBs are spinning at periods that are longer 
than expected at equilibrium, unless the magnetic fields 
are stronger than 10^^ G (c.f., Santangelo et al. 2011). 

We continue to observe an asymmetry in the distri- 
bution of absorbed HMXBs in the inner spiral arms 
(Bodaghee et al. 2007). Figure 3 presents HMXBs in 
the central 150° of the Milky Way where the symbol 
size reflects the reported A^h value. Within 30° of the 
Galactic Genter, there is a clear right-left asymmetry 
with 4 times as many absorbed sources at negative lon- 
gitudes (right half of the figure) as there are at positive 
longitudes (left half of the figure). Absorbed sources 
are clustered in the direction of the Inner Perseus and 
Norma Arms with no equivalent overdensity towards the 
other arms. The Norma Arm happens to be the direction 
which features the highest number of massive-star form- 
ing regions (Russeil 2003). These are enormous clouds 
whose fragments spawn the massive precursor stars that 




Fig. 4. Galactic distribution of HIVIXBs whose distances are l<nown (75, filled triangles) and the locations of OB star-forming complexes (464, 
circles) from Russeil (2003). The symbol size of the latter is proportional to the amount of activity in the complex. The spiral arm model 
of Vallee (2008) is overlaid with the Sun situated at 7.6 kpc from the center. HMXBs whose distances are not known have been placed 
at 7.6 kpc (26, empty triangles). The shaded histogram represents the number of HMXBs in each 15° bin of galactic longitude. 



eventually become HMXBs. In the next section, we will 
examine their spatial link in detail. 

3. Spatial Correlation 

One of the goals of statistical analyses of HMXB popu- 
lations is to compare their spatial distribution with that 
of the massive-star forming regions in which they are 
suspected of being born. Thus far, statistical compati- 
bility between HMXBs and young star-forming systems 
has been demonstrated in distributions of galactic longi- 
tudes or galactocentric distance. The drawback is that 
these histograms consider only a single dimension. An 
entire axis is projected onto the other entailing infor- 
mation loss. This obstacle can be overcome by consid- 
ering the spatial (or cross-) correlation function i^. For 



a given HMXB in a volume element SVi, the probabil- 
ity SP (above Poisson) of finding an OB complex in a 
volume element SV2 separated by a distance r is: 

SP = nin2[l+^{r)]5ViSV2 (1) 

The mean number densities are given by rii and 
712 where the subscripts refer to the two populations: 
HMXBs (1) and OB complexes (2). Young objects such 
as HMXBs and OB stars are formed near the disk, i.e., 
we can safely assume that the displacement along the 
z-axis is negligible, so surface (rather than volume) ele- 
ments are used. Around each HMXB, we construct rings 
of 1 kpc thickness. In each ring, we count the number of 
observed OB complexes, and we keep a separate tally of 
the number of OB complexes from a random distribu- 
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Fig. 5. Spatial correlation functions ^(r) from both estimators are presented for three randomization models: a Gaussian ring centered at 
7.6 kpc with (Tz = 2; a uniform disk for < r < 16 kpc; and a uniform disk with a hole for 2 < r < 16 kpc. Error bars represent 3-ct 
uncertainties. At short distances from a given HMXB, ^ is greater than 1. This indicates that the closest neighbors of an HMXB tend to 
be observed OB associations rather than OB associations drawn from a random distribution. 



tion. An HMXB-OB pair is referred to as a I?L'i2-pair 
(for data-data pair: the numbers designate the popula- 
tions), while an HMXB-random OB pair is labelled DR12 
(for data-random). In this way, we construct ^ for each 
radius according to the definition of Peebles (1980): 



estimator for ^ in which the variance is nearly Poissonian: 



nRDDi2 

4(r) — 1 

nDDRi2 



(2) 



Clearly, if ^ = 0, which implies that each ring contains 
as many DD-pairs as DR-pairs, then Eq. 1 is simply a 
uniform Poissonian probability. However, if ^ > 0, then 
there is a higher chance of an HMXB having a neigh- 
bor that is an observed (rather than a randomized) OB 
complex. Landy & Szalay (1993) propose a more reliable 



DD12 — DR12 — DR21 + RR12 
RR12 



(3) 



Figure 5 shows that the neighbor of a randomly- 
selected HMXB is more likely to be a known OB complex 
than expected from Poisson statistics. For r < 3 kpc, 
the clustering signal from the Peebles (1980) estimator 
is 7-llf7 in excess of Poisson (up to 17(t for the Landy 
& Szalay (1993) estimator). The observed and random 
surface-density distributions are statistically compatible 
only at large radii (i.e. r > 5 kpc). 

In other words, HMXBs and OB associations are clus- 
tered together, as one would expect since HMXBs are 
young systems that have not ventured far from their 



likely birthplaces. This is the first time the relation- 
ship has been demonstrated in Cartesian space in our 
Galaxy. The characteristic scale of ^ holds clues to the 
average offset of the HMXB population compared to the 
current crop of star-forming regions. In turn, this will 
help constrain the amount of migration and kinematical 
timescales of HMXBs. 

It is reasonable to expect symmetry in the distribution 
of HMXBs and OB complexes around the galactic disk 
and our randomization models assumed this. However, 
the fact that there are less HMXBs detected (with dis- 
tances measured) behind the Bulge will surely add bias 
to ^. Therefore, we generated ^ for a restricted bound- 
ary corresponding to a circle of radius 8 kpc around 
the Sun — a perimeter inside of which most HMXBs and 
OB complexes should be detectable and their distances 
known with reasonable accuracy. When only pair counts 
of objects (53 HMXBs and 363 OB complexes) in the So- 
lar neighborhood (r < 8 kpc) were considered, the clus- 
tering signal at small radii persisted with a statistical 
significance of lo for r < 1 kpc. 

Clustering between HMXBs and OB complexes is ex- 
pected. How would ^ react if HMXBs were compared 
to a population of sources for which no spatial corre- 
lation is anticipated? We tested the correlation func- 
tions of HMXBs against a set of 133 globular clusters 
from Bica et al. (2006) that are located less than 20 kpc 
from the GC. Globular clusters contain older populations 
such as cool KM dwarf stars and low-mass X-ray bina- 
ries (LMXBs). Unlike HMXBs, globular clusters (and 
LMXBs) are densely packed in the Galactic Bulge and 
their numbers drop exponentially with increasing radius 
from the GC. The random distribution of globular clus- 
ters was thus modeled as an exponential decay law ad- 
justed to be consistent with the observed distribution. 
The ^ between HMXBs and globular clusters is consis- 
tent with for r < 1 kpc from an HMXB. In other words, 
for any given HMXB, its immediate neighbors were just 
as likely to be drawn from the observed set as from the 
random set. Shot noise dominates at large distances. 

Currently, we are working on accounting for other ob- 
servational biases, and we are examining the behavior of 
^ under different evolutionary timescales. Further dis- 
cussions of these results will soon be submitted for pub- 
lication as Bodaghee et al. (2011b). 



confirm theoretical predictions (such as the link between 
HMXBs and OB associations), and to uncover trends 
that await explanation (such as the correlation in the 
A^H vs. spin period diagram). With large and growing 
samples representing each of the three known HMXB 
classes, we can achieve better statistics while exploring 
the parameter space occupied by HMXBs. A synthetic 
study of the different populations of HMXB will help 
reveal clues to the nature and evolutionary history of 
each class, which will help position each group within a 
unified model of HMXBs. 
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4. Conclusion 

Multi-wavelength analysis of an HMXB can be used to 
reveal the mechanisms responsible for the emission, de- 
scribe the composition and morphology of material sur- 
rounding (and in some cases enshrouding) the system, 
and permit a firm identification of the optical/IR coun- 
terpart. Multi-variate analysis of the spatial, timing, and 
spectral parameters of HMXB populations has served to 



